Exotic singularities and spatially curved Loop Quantum Cosmology 
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We investigate the occurrence of various exotic spacelike singularities in the past and the future 
evolution of k — ±1 Friedmann- Robertson-Walker model and loop quantum cosmology using a 
sufficiently general phenomenological model for the equation of state. We highlight the non-trivial 
role played by the intrinsic curvature for these singularities and the new physics which emerges at 
the Planck scale. We show that quantum gravity effects generically resolve all strong curvature 
singularities including big rip and big freeze singularities. The weak singularities, which include 
sudden and big brake singularities are ignored by quantum gravity when spatial curvature is negative, 
as was previously found for the spatially flat model. Interestingly, for the spatially closed model 
there exist cases where weak singularities may be resolved when they occur in the past evolution. 
The spatially closed model exhibits another novel feature. For a particular class of equation of 
state, this model also exhibits an additional physical branch in loop quantum cosmology, a baby 
universe separated from the parent branch. Our analysis generalizes previous results obtained on 
the resolution of strong curvature singularities in flat models to isotropic spacetimes with non-zero 
spatial curvature. 

PACS numbers: 04.60.Pp,04.20.Dw,04.60.Kz 



I. INTRODUCTION 



A fundamental limitation of the cosmological models 
based on general relativity (GR) is the occurrence of sin- 
gularities. Perhaps one of the simplest examples is the 
case of an expanding homogeneous and isotropic universe 
filled with a matter satisfying strong energy condition 
such as dust or radiation. Independent of the intrinsic 
geometry of the universe, be it closed, flat or open, the 
past evolution of such a universe from arbitrary initial 
conditions leads to an initial singularity: the big bang, 
where the classical dynamical equations break down and 
the physics stops. Another example is the case of in- 
flationary universe in which even though the evolution 
is almost deSitter, the classical spacetime is past incom- 
plete 

In recent years, various new singularities have been 
found in classical cosmology Unlike the big bang 

and big crunch singularities, these singularities do not 
occur when scale factor vanishes. These occur either at 
finite values of the scale factor or when it diverges. Recall 
that for the matter satisfying weak energy condition, the 
latter is not possible as the spacetime curvature goes to 
zero when scale factor goes to infinity. However, if mat- 
ter violates weak energy condition, for example in the 
case of phantom fields, then spacetime curvature will di- 
verge as scale factor becomes infinite. For homogeneous 
and isotropic models with matter equation of state in the 
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form of a perfect fluid these exotic singularities come in 
four types. Big rip (type I) where the energy density 
and pressure diverge along with a divergence in the scale 
factor, sudden singularity (type II) occurring at a finite 
value of the scale factor and energy density with a di- 
vergence in pressure, big freeze (type III) where energy 
density and pressure diverge at a finite value of the scale 
factor and big brake (type IV) singularity where scale 
factor, energy density and pressure are finite but there 
is a divergence in the time derivative of the pressure or 
rate of change of energy density. 

Lack of successful resolution of these singularites ren- 
ders classical cosmological models incomplete. It is gen- 
erally believed that existence of these singularities is a re- 
sult of assuming the validity of GR even in the regime of 
large spacetime curvature where the effects due to quan- 
tum gravity may become important and lead to signifi- 
cant departures from the classical theory. It is thus hoped 
that incorporation of quantum gravitational effects may 
result in a possible resolution of these singularities. 

Loop quantum gravity (LQG) is one of the candidate 
theories of quantum gravity which attempts to address 
this issue. It is a non-perturbative and background in- 
dependent quantization, with a key prediction that the 
continuum differential geometry of the classical theory is 
replaced by a discrete quantum geometry in the quan- 
tum theory. Perhaps one of the best illustrations of the 
novel effects of quantum geometry is captured in loop 
quantum cosmology (LQC) which is a quantization of 
homogeneous spacetimes based on LQG A key 

prediction of LQC is that the big bang singularity is re- 
placed by a big bounce, which is a direct consequence of 
the underlying quantum geometry [§-12]. These results 



which were first obtained for homogeneous and isotropic 
models (for all values of spatial curvature) with a mass- 
less scalar field have been extended to inflationary po- 
tential [lU , anisotropic spacetimes [3 UH and also cer- 
tain inhomogeneous situations [lq ]. Further, using an 
exactly solvable model it has been shown that the ex- 
pectation values of energy density have a universal up- 
per bound for a dense subspace in the physical Hilbert 
space [17[ • There are strong constraints on the change in 
relative fluctuations of quantum observables across the 
bounce [Isj . Recently, much stronger constraints on the 
change in dispersions have been obtained by Kaminski 
and Pawlowski These results show that a universe 
like ours i.e. macroscopic at late times bounces from a 
a similar universe at very early times (in the contracting 
branch) and the universe recalls almost of all its state 
through the bounce. 

Interestingly, the loop quantum dynamics admits an ef- 
fective description on a continuum spacetime which can 
be obtained using coherent state techniques [2(| Hl[ . An 
important feature of this analysis is that one can ob- 
tain an effective Hamiltonian from which one can obtain 
modified Fricdmann and Raychaudhuri equations as the 
Hamilton's equations. The modified set of dynamical 
equations inherit quantum geometric effects via higher 
order non-perturbative corrections which vanish at small 
spacetime curvatures. It is important to note that vari- 
ous numerical simulations have shown that effective equa- 
tions capture the underlying quantum evolution very ac- 
curately for universes which become macroscopic at late 
times. These thus prove to be useful tools to understand 
the physics in LQC, such as whether the underlying the- 
ory has well defined ultra-violet and infra-red limits. It 
turns out that even though there exist various quantiza- 
tion ambiguities, there is a unique quantization lea ding 
to a consistent unambiguous physical description |22l . |23[ 
(the improved dynamics [Til Il2j: which is being consid- 
ered here). 

Using effective equations, we can ask various questions 
regarding the generality of singularity resolution in LQC. 
For example, one can ask whether spacetime curvature 
is always bounded in LQC? Here we should note that 
a universal bound on energy density (as in LQC), does 
not imply that the spacetime curvature is also bounded. 
This is easy to understand for the classical cosmology, 
where the Ricci scalar, which provides us a complete in- 
formation about the spacetime curvature in the homo- 
geneous and isotropic spacetime, depends both on the 
energy density and pressure. Though for most matter- 
energy configurations, the behavior of equation of state is 
such that an upper bound in energy density is sufficient 
to control the divergence in pressure and hence the space- 
time curvature, it is not difficult to come up with counter 
examples with a more general equation of state [U [Hj| ■ 
Hence, an upper bound in energy density is not sufficient 
to prevent a divergence in the spacetime curvature. 

A pertinent question is whether this divergence sig- 
nals the end of spacetime in LQC. In order to answer 



this question, we recall that even in GR we encounter 
events where spacetime curvature blows up but there is 
no associated physical singularity. This can happen if 
the tidal forces are not strong enough to cause a com- 
plete destruction of in-falling objects in to the singular- 
ity and geodesies can be extended beyond such events. 
It turns out that the events where spacetime curvature 
diverges in flat isotropic LQC are weak singularities and 
geodesies can be extended beyond them. In flat isotropic 
LQC, the divergence of spacetime curvature occurs only 
for sudden singularities which are caused by a divergence 
in pressure at a finite scale factor and energy density. It 
is straight forward to show that the expansion parameter 
in this case is bounded and the spacetime is geodesically 
complete in the flat isotropic LQC pBj . 

In this paper we take the first step to generalize the 
above result by including intrinsic curvature in the space- 
time. This is done by considering the effective dynamics 
of loop quantized spatially closed and open models in 
the Robertson- Walker geometry. In the classical Fried- 
mann dynamics, intrinsic curvature term enters in form 
of 1/a 2 term in the dynamical equations. Thus it is 
quite straightforward to understand the expected mod- 
ifications from the results in the spatially flat model. 
In LQC, the quantization of intrinsic curvature brings 
non-trivial modifications to the effective description and 
makes the resulting form of effective dynamical equations 
less straight forward to analyze. Though one expects that 
at small intrinsic curvatures one recovers the results of 
flat isotropic LQC, physics may be bring up surprises 
when intrinsic curvature is large. As we will show, this 
is indeed what happens in the case of k = ±1 models in 
LQC. 

Our analysis will be based on considering a sufficiently 
general phenomenological model for the equation of state 
which was proposed in Ref.[26|. This model allows a 
study of all of the exotic singularities by the choice of 
different parameters and was earlier used for investiga- 
tion of resolution of strong curvature singularities in the 
flat isotropic LQC 25]. As we will see, the effective 
dynamical equations of spatially curved models approxi- 
mate those of the flat model in LQC at large volumes be- 
cause the contribution from intrinsic curvature becomes 
negligible in this limit. Thus for future singularities, at 
large volumes, the resulting physics is similar for models 
with or without spatial curvature in LQC. However for 
certain values of parameters the spatially closed model 
in LQC permits two separate physical branches, a short 
lived baby universe at small volume and a parent uni- 
verse which evolves to a macroscopic size. This branch 
is absent in the classical theory and has a pure quantum 
geometric origin. 

To completely capture the new physics from inclusion 
of intrinsic curvature, it is important to study exotic sin- 
gularities in the past evolution when they occur at small 
volumes. Our analysis of past and future exotic singu- 
larities shows that all strong curvature singularities are 
resolved in k = ±1 isotropic LQC. The scale factors at 



2 



which big rip and big freeze singularities occur are ex- 
cluded from the allowed range by loop quantum effects. 
As in the flat model, the spacetime curvature can diverge 
in spatially curved LQC, however when ever this happens 
one has a weak singularity which is known to be harm- 
less. In almost all cases these singularities are ignored 
by LQC. The only exception to this occurs in spatially 
closed model where weak singularities occurring in the 
past evolution may be resolved. This occurs purely be- 
cause of the non-trivial role of intrinsic curvature effects 
in LQC. 

We organize our paper as follows: In Sec|Tl] we re- 
visit the classical equations for spatially curved model 
in classical cosmology and introduce the phenomenolog- 
ical ansatz of the equation of state. To facilitate the 
reader to follow the derivation of effective equations in 
LQC, we derive the Friedmann and Raychaudhuri equa- 
tions in the Hamiltonian framework. In Sec lIIIl we derive 
the effective equations for spatially curved LQC starting 
from the effective Hamiltonian [lj, H3, [H| • Using these 
modified equations, we numerically obtain solutions and 
discuss the new physics in Sec lIVl Here we show that 
all exotic strong curvature singularities, irrespective of 
whether they occur in the past or the future, are resolved 
in spatially curved LQC. We summarize our results in 
SecIV] 



II. CLASSICAL THEORY 
A. Hamiltonian cosmology 

The fundamental equations of cosmology describe the 
evolution of the scale factor a in proper time t. These 
equations can be derived in a simple Hamiltonian frame- 
work. To this purpose, consider the conjugate variables 
(c, p) where c is the Ashtekar-Barbero connection and p 
is the triad (which without any loss of generality will be 
chosen with positive orientation). These are related to 
the metric variables as 



c = 7a + k 



P : 



(1) 



with the relation between c and a valid only in the classi- 
cal theory and it modifies when quantum gravity effects 
are considered. Here d = da/dt, k = 0, ±1 is the normal- 
ized intrinsic curvature and 7 £ R is the Barbero-Immirzi 
parameter 1 . The conjugate variables satisfy the following 
Poisson bracket 



In these variables, the Hamiltonian for gravity in a ho- 
mogeneous and isotropic spacetime becomes 



Ug 8ttG 



7 ~ 2 ^ [(c-fc) 2 + fc 7 2 



(3) 



We introduce a generic matter field with Hamiltonian 
T~im = P 5 Pi where p is the matter energy density. The 
total Hamiltonian of the system % = H g + T-L m is con- 
strained to vanish. 

The dynamics is given by the Hamilton's equation 

p = { Pi n} = -{c,p}— = 2 7 - 1 ^p(c-k) . (4) 

Since the Hamiltonian is a constraint, namely % ~ 0, we 
find that 

(c- fc) 2 = ^^-~i 2 pp - fc7 2 . (5) 
We can now recover the usual Friedmann equation as 



H 



(6) 



In order to obtain the Raychaudhuri equation for the 
acceleration, we compute the equation of motion for c 



c = {c,H} = {c,p} 



dU 
dp 



+ 



(c - fc) 2 + fc 7 2 ( 8irG-fdH 
47rG7 



3 dp 

VP(P + 3P) (7) 



where we have introduced the thermodynamic pressure 
P as the derivative of Ti. m with respect to the volume. 
This gives dH m /dp = —\^fp P- From (J4j) and dU) we 
compute 



.. IGirG oI 2 _ . 
p = ——pp-2k + -ypp c 

3 7 



and we obtain 2 



H = = ^ (p - 3P) 



2p ■' [, 
that yields the Raychaudhuri equation 
a 



H + H 



\irG(p + 3P) 



(8) 



(9) 



(10) 



. , 8ttG 
i c , Pi = ^—7 



(2) 



1 The value of 7 can be fixed by computing the black hole entropy 
in LQG. In our numerical analysis we set 7 fs 0.2375 1291 . 



2 Equivalently, one can obtain H using 2HH = ^f^-p + 2H^ and 
p = —3H(p + P). The conservation equation will hold also in 
the case with quantum gravity modifications, since it results from 
the Hamilton's equation for the matter part, without involving 
the gravitational part. 
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Note that a divergence in p and/or P can lead to a di- 
vergence in the Ricci scalar 



R = 6 (H + 2H 2 



(11) 



Analogously, a divergence in P, the derivative of the pres- 
sure, can lead to a divergence in the derivative of the 
curvature 



R = 6 ( H + 4HH - 2H- 
P 



(12) 



Divergences can similarly be computed for the higher 
derivatives. 



B. Phenomenological model 



singularity 


divergences 


parameters 


Type I 


a — > oo, p — ¥ oo, P — > oo 


3/4 < a < 1 Wl, B 


Type II 


P — > oo 


a < or A/B > 1 


Type III 


p — > oo , P — > oo 


a > 1 \M, B 


Type IV 


P — > oo 


0<a<l/2 VA,B 



the value of A and P. No exotic singularity appears for 
1/2 < a < 3/4. 

In Sec lIVI we study these singularities, in particular 
by showing the behavior of the Hubble rate, the Ricci 
curvature and its derivative, since these quantities are 
closely related to the energy density, the pressure and its 
derivative. 



Different types of singularities are categorized depend- 
ing upon whether the divergences appear in the scale fac- 
tor a, in the energy density p, in the pressure P or in its 
derivative P. For our investigation we choose a specific 
expression for the pressure 



P 



-P- 



ABp 



2a- 1 



Ap° 



B 



(13) 



that allows to obtain various singularities by varying the 
parameters of the model A, B and a (26|. The derivative 
of the pressure is given by the expression 



P = 3H(p+P) 



1- 



(2a— 1) AB p 2a 
Ap a - X + B 



(I- a) A 2 Bp' 



3q-3 



VBf 
(14) 

The equations of motion of the model can be integrated 
and the scale factor can be expressed as a function of the 
energy density p as 



a = a a exp 



{2A + Bp 2 )p 
6AB(l-a) 



(15) 



Here a is an integration constant, and for the singu- 
larities at finite scale factor, it corresponds to the value 
where the singularity appears. Inverting this expression 
we obtain 




V(i-") 



(16) 



The following table summarize the relation between the 
parameters of the model and the quantities that diverge, 
while the other quantities remain finite. 

In the classical theory, singularities can appear in the 
past or in the future, depending on the choice of the 
parameters A and B. If A and B have different sign, 
they always give rise to a sudden singularity. The other 
singularities depends on the value of a, irrespectively on 



III. EFFECTIVE DYNAMICS IN LQC 

The underlying dynamics in LQC captures quantum 
discreteness of LQG. However, for coherent states one 
can obtain an effective continuum spacetime description. 
Such an analysis has been carried for different types of 
matter and an effective Hamiltonian has been obtained 
[20l l2~j| . The effective Hamiltonian provides effective dy- 
namics via Hamilton's equations up to controlled higher 
order approximations. Extensive numerical simulations 
show that the effective dynamics is an excellent approx- 
imation to full quantum dynamics for states which cor- 
respond to a macroscopic universe at late times [lll - fl3l |. 
The non-local quantum gravitational effects originating 
from the underlying quantum geometry primarily modify 
the gravitational part of the Hamiltonian. The matter 
part of the Hamiltonian constraint remains unaffected. 
(We discuss the underlying approximation later in a re- 
mark) . 

The effective Hamiltonian can be written as [TTI . [i~2l . 
Sill 11 11 : 



off ■= 



A(v) 
16ttG 



where 



p p 



(p{c-k)) - k X 



M 



(17) 



(18) 



and x nas different expressions for k = 1 and k = — 1 
geometries. For k = 1 it is given by [l2[ 



X := sin 2 p, - (1 + j 2 )p 2 
and for k = -1 it becomes H 11 



2 -2 

X ■= -7 M 



(19) 



(20) 



In above equations, A denotes the minimum eigenvalue 
of the area operator in LQG and v denotes the eigenval- 
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ues of the volume operator V = p 3 ^ 2 in LQC. 3 Without 
any loss of generalization we restrict ourselves to the pos- 
itive eigenvalues of the volume oeprator. 



(21) 



with K = 2/3v / 3V3. 

For v > 1, the expression for A(v) yields [TTJ] 



A(v) 



27K£ P1 [8tc. 



2 7 3/ 2 



6 



dp 1 ' 2 
p 2 7 2 



(22) 



where in the last step we have used Eq.([2T|). Thus, the 
effective Hamiltonian becomes 

H cS = -g^2 (sin 2 (p(c - k)) - k X ) + Hu ■ (23) 

The vanishing of the Hamiltonian constraint, 1-L c ff ~ 0, 
leads to 



sin 2 (^(c — k)) = 



8ttG YP- 



2 r,2 



M 



kx 



k X (24) 



3 y/p " Pcrit 

where we have defined the critical energy density p cr it as 



Pcrit 



87rG7 2 A 



(25) 



The modified Friedmann equation can be obtained 
from the effective Hamiltonian (|23l) . by computing the 
Hamilton's equation for p: 

2 _ 

p = {p,H cS } = — s/p sin (p(c-fc)) cos (p(c-fc)) . (26) 
The Hubble rate H 2 = a 2 /a 2 , then becomes 



H 2 



2p J 7 2 A 



sm 



(/2(c-fc))(l-sin 2 (p(c-fc))) 



8 ~ , fc X \ / p 
-vrGp + -5— 1 fex 

3 7 A/ V Pcrit 



This equation can be rewritten as 



(27) 



F 2 = ^(p-Pr)(— (p 2 -p)) (28) 

3 V Pcrit 



where p\ and p 2 are defined as 

. = ~3fcx 
Pl ' 8ttG 7 2 A " 



(29) 



3 Note that the value of K is different from various previous works. 
This is due to change in the expression for the area gap in LQC 
1 1 5| . In comparison to Ref. [ill , we have p. 2 p 2 = &i/3wy£p\. 



and 

P2 ■= Pcrit (1 - kx) ■ (30) 
In the classical limit, A — > 0, we obtain 



X 



3k 1 
-7 2 p 2 , pi -> and - — (p 2 -p) -> 1 . (31) 



Pcrit 



Thus we recover back the classical Friedmann equation 
in the limit A — >• 0. 

Similarly, using the Hamilton's equation for c, we can 
obtain the modified Raychaudhuri equation: 

2 47rG / o 167rG f P , V 3 „ 

a 3 3 VPcrit / \ 2 

+ J ?a(— + ^)-^(— + fc X-|>2) 

7 2 A V Pcrit / 7 A V Pcrit 2 7 



where 



£ := sin 2 (/i) — /i sin(p) cos(p) 



(33) 



Using eq. ([32l) and eq. ([2T|) we obtain the equation for the 
rate of change of the Hubble rate 

H=(-4,G( P + P) + k -^-) (l-2(j^ + kx 



7 2 A 



Pcrit 



(34) 

It is then straightforward to verify that the conservation 
law 



p + 3H(p + P) =0 



(35) 



remains unchanged in the effective dynamics of LQC. 

We can now write down the expression for the Ricci 
scalar which captures the complete behavior of the space- 
time curvature for the homogeneous and isotropic model: 



R 



6 H 2 



A-nG , , 87rG ( p , . 

^Hp-3P) + ^t- \— + k X ){p + 3P) 



3 V Pcrit 
P 



k^_ 2Cfc (_p_ +ky _l\ k 
7 2 A 7 2 A \ Pclit + X 2 J + a 2 



(36) 



Further, as for the case of modified Friedmann equa- 
tion, the modified Raychaudhri (|3"2"1) . H ([M]) equations 
and Ricci scalar go to the classical GR versions in the 
limit A -> 0. 

It is important to point out a notable difference be- 
tween the Ricci scalar in the classical theory and LQC. 
In the classical theory, the expression for R is indepen- 
dent of the curvature index k 



R = 8ttG(p - 3P) 



(37) 



Whereas in LQC, Ricci scalar is dependent on the 
curvature index. The expression for the Ricci scalar is 
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different for flat, open and closed models in LQC. We 
will see in the next section that it leads to interesting 
distinctions for singularity resolution for different spatial 
geometries in LQC. 

We conclude this section with the following remark: 

Remark: In the derivation of effective equations we have 
worked under the approximation that contributions from 
the inverse scale factor effects are negligible. The approx- 
imation is well motivated due to two reasons. Firstly, 
from insights gained from various numerical simulations 
which show that in LQC, modifications originating from 
the non-local field strength of Ashtekar-Barbero con- 
nection, which results in the trigonometric function in 
eq. (|17[) . significantly overwhelm the modifications origi- 
nating from inverse scale factor [III [13, H3. Secondly, 
for v > 1 inverse scale factor effects are in any case negli- 
gible (in the fundamental representation of the theory ). 
This corresponds to a > l.5£p\ from (|2"Tjl . In our analysis 
we will only consider scale factors large than this value. 
Thus in our analysis it will be safe to make this approxi- 
mation. We expect the results to hold true in general as 
inverse scale factor effects tend to weaken the strength 
of gravity (or make it effectively more "repulsive" ) , thus 
aiding singularity resolution. 



IV. PHENOMENOLOGICAL MODEL: 
NUMERICAL RESULTS 

The phenomenological model introduced in Sec. IIB 
allows us to study different kinds of singularities by 
choosing appropriate values for various parameters in 
Eq. (|13[) . If A and B have a different sign, then there is 
always a sudden singularity independently of the value 
of a. If A and B have the same sign, then we obtain the 
following classification: 

Type I singularity (Big Rip): A > and 3/4 < a < 1. 

Type II singularity (Sudden): a < 

Type III singularity (Big Freeze): a > 1 

Type IV singularity (Big Brake): < a < 1/2 

If a — or 3/4 < a < 1, there are no type I-IV 
singularities. Further, apart from the type I singularity 
which occurs with an associated divergence in the scale 
factor, type II-IV singularities can occur both in the 
past or the future of an expanding branch in the classical 
FRW model. 



4 The situation changes when one is considering higher represen- 
tations where inverse scale factor effects can lead to novel phe- 
nomenological effects (see for eg. R.cfs. 30 32]). See, however 
Ref. [33ll for subtleties in dealing with higher representations. 



Using the effective loop quantum dynamics as elabo- 
rated in the previous section, we carried out extensive 
numerical simulations for various parameters with the 
phenomenological equation of state (flB"]). Below we dis- 
cuss various exotic singularities and show the represen- 
tative results for different cases. (We use Planck units to 
show results in various plots). 

1. Type I singularity: The Big Rip 

A type I singularity is also called "Big Rip" because 
there exists a finite time in which the scale factor, energy 
density and pressure diverge, tearing apart the universe. 
The dominant energy condition is broken and the equa- 
tion of state converges to w — — 1 when approaching the 
singularity. 
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FIG. 1: Type I future singularity: Hubble rates for 
k = 0,±1 models are shown. Dashed curves corresponds to 
the classical theory and the solid ones to LQC. (Curves for 
different values of curvature index coincide with each other for 
the resolution in the plot). In the classical theory, the Hubble 
rate diverges in a finite time with a —¥ oo. We see that the 
effective LQC curves behave similarly to the classical ones for 
small values of the scale factor a when the energy density is 
very small compared to the Planck scale. As the energy den- 
sity increases, departures from classical trajectories became 
significant and we see a quantum recollapse. In this figure 
the parameters are a = 1000, A = 0.1, B = 1 and a = 0.8. 

In LQC, the energy density is bounded by a maximal 
value, therefore all type I singularities are resolved by the 
quantum theory. The energy density grows as in the clas- 
sical theory as the singularity is approached, but when it 
reaches close to Pent, quantum effects lead to significant 
modifications to the dynamical trajectory. The accelera- 
tion a becomes negative and the Hubble rate goes to zero 
(Fig.[T]). Instead of ripping apart in finite time, the loop 
quantum universe recollapses and the evolution contin- 
ues. The presence of a maximal density affects also the 
curvature invariants. In particular the Ricci scalar and 
its derivatives remain bounded during the whole evolu- 
tion. 

This behavior which was first noted in k — model (2f| 
(see also Ref. [34| IH[ for earlier results) remains essen- 
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FIG. 2: Zoom of the Hubble rate around a„ = 1 000. The 
curves are respectively thin (blue) for k = — 1 on the left, 
black for k = and thick (red) for k — +1 on the right. 
The parameters are the same as in Fig.[l] This zoom reveals 
the differences for evolution in Hubble rates in the region 
around a , that corresponds in this case to the "initial point" 
of cosmic evolution. 



tially the same in the k = 1 and k = — 1, as is shown by 
Fig.Q] This is not surprising because type I singularities 
occur in future in the expanding branch where the scale 
factor is very large and the effects due to intrinsic cur- 
vature are expected to be negligible. (We will see below 
that effects due to intrinsic curvature can indeed produce 
surprising results for other exotic singularities). 

The only notable difference between the k = and 
k = ±1 cases is that the exact value of the density at 
the bounce changes, but the correction is small (of the 
order of ^pi/a^ in ) and has no effect on the qualitative 
behavior of the trajectory. The three cases differ only 
around the value a of the scale factor as shown in Fig. [5] 
As one would expect for this equation of state, at large 
scale factors the evolution in k = ±1 model mimics that 
in the flat model and the type I singularity is resolved. 
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FIG. 3: Type II future singularity: Hubble rate and 
Ricci scalar for k — in black, k — + 1 in thick (red) and 
k = — 1 in thin (blue) curves are shown. Classical curves 
are dashed and LQC curves are solid. The Hubble rate goes 
to zero both in the classical theory and LQC for all values 
of the curvature index at a = a = 100. The Ricci scalar 
diverges at a in the classical theory and LQC. (Since Ricci 
scalar is independent of the value of curvature index in the 
classical model, there is only curve for the classical theory). 
For k = +1 there exist two disjoint solutions of the Friedmann 
equation, and the further branch appearing there is bounded 
both in the Hubble rate and in the spacetime curvature. The 
parameters are A = — 0.1 , B = 10 and a = 1/4. 



2. Type II singularity: The sudden singularity 

In the type II singularity, the energy density does not 
diverge. The singularity is caused by a divergence in 
the pressure which results in a divergence in spacetime 
curvature at a finite value of the scale factor. To under- 
stand this singularity in more detail, let us consider Eq. 
(1131) : the pressure diverges when the energy density is 

p = p s = (-A/B) . Inserting this expression into 
(fl~5|) we find that the singularity appears for a — > a a . 
It is important to note that in LQC, quantum geomet- 
ric effects do not regulate any divergence in pressure and 
the spacetime curvature can diverge p5| . It turns out 
that such a divergence for the type II singularity does 
not signal the end of spacetime. Geodesies can be ex- 
tended beyond this singularity. Sudden singularity also 
turns out to be a weak singularity because the tidal forces 
are not strong enou gh t o cause a complete destruction of 
arbitrary detectors 25]. Thus weak singularities signal 



neither a break down of the underlying theory nor the 
end of spacetime. 

In Fig. |3l we show the evolution of the Hubble rate for 
different values of spatial curvature index for type II sin- 
gularity occurring in the future of an expanding branch. 
We see that both in the classical theory (dashed curves) 
and LQC (solid curves), the Hubble rate vanishes in a 
finite time for all values of the spatial curvature at scale 
factor a — a . In contrast, the plot of the Ricci scalar 
shows that it diverges at this value of the scale factor for 
all values of the spatial curvature in the classical theory 
as well as in LQC. 

As is clear from Fig.[3j the results about future sudden 
singularity are on expected lines of the k = model when 
the scale factor is very large, where the effects due to in- 
trinsic curvature become negligible. However, a surpris- 
ing result emerges for the k = 1 model. Our numerical 
analysis shows that, for certain values of the parameters, 
there exists an additional baby evolutionary branch sep- 
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FIG. 4: Type II past singularity: Hubble rate and 
Ricci scalar for k — and k = ±1 are shown. Classical 
curves are dashed and LQC curves are solid, k = is in 
black, k — +1 is thick (red) curve and k — —1 in thin 
(blue) curve. The singularity is at a — 2, where the 
Ricci scalar diverges. Note that for k = +1 the universe 
begins at a > a , and the Ricci scalar is always finite. In 
this figure the parameters are A = 0.1, B — —10 and a = 1/4. 
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FIG. 5: Type III past singularity: Classical (dashed) and 
effective LQC (solid) Hubble rate and Ricci scalar are shown. 
The curves are respectively thick (red) for k = +1 starting 
on the left, black for k = in the middle and thin (blue) for 
k = — 1 stating a bit more on the right. In the classical case 
there is a divergence of both H and R at a — 1 while the 
curves for LQC remain bounded (in the figure the Ricci scalar 
for k = — 1 is covered by the black line of the k — case). 
Here the parameters used are A = —100, B = — 1 and a — 2. 



arated from the main evolutionary branch. For k = +1, 
H 2 is positive in two disjoint intervals. Consequently 
there exist two disjoint branches corresponding to real 
solutions of the Friedmann equation. The extra branch 
occurs at small scale factors and signifies the non-trivial 
new physics which emerges from the quantization of in- 
trinsic curvature in LQC. The new branch, that we can 
see in the Fig.(3J is bounded both in the Hubble rate and 
in the spacetime curvature. 

Another interesting feature which can be seen from 
Fig. [3] is that for k = 1 model, the classical curve (shown 
with a dotted curve) depicts that there is no initial sin- 
gularity. The Hubble rate vanishes both in the past and 
the future of the classical evolution. The classical uni- 
verse faces a sudden singularity in future, but since it is 
a weak singularity, the classical k = 1 universe is geodesi- 
cally complete for the considered equation of state. This 
does not hold true for the k = and k — — 1 case, and 
these classical spacetimes are past incomplete. 

Fig.|4]shows the plot of Hubble rate and Ricci scalar for 
different values of curvature index when the sudden sin- 
gularity occurs in past. In such a case the effects due to 
intrinsic curvature lead to novel features in the physics of 



singularity resolution. We find that classically the Hub- 
ble rate vanishes for all values of k at a — a and the 
Ricci scalar diverges. However in LQC, though the Hub- 
ble rate vanishes for values of k, the divergence in Ricci 
scalar does not occur for a closed universe. It turns out 
that as scale factor approaches a , the Ricci scalar in- 
creases to large values but it does not diverge for the al- 
lowed values of the scale factor. (The divergence occurs 
in the forbidden region of effective LQC dynamics where 
the Hubble rate is imaginary). This result implies that 
in the closed model for the parameters considered here, 
quantization of intrinsic curvature leads to a resolution 
of the past sudden singularity. This is a direct conse- 
quence of the way intrinsic curvature terms enter in to 
the expression of the Ricci scalar in LQC. We emphasize 
that type II singularities whether in future or past, even 
though they are harmless, are not resolved in LQC for 
k = and k = — 1 models. In contrast we find that these 
singularities are resolved, for a certain choice of parame- 
ters, when occurring in the past for k = 1 model. 
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3. Type III singularity: The Big Freeze 

A type III singularity occurs at a finite value of scale 
factor where both the Hubble rate and the Ricci scalar di- 
verge. In LQC, the presence of a maximal energy density 
Pcrit bounds the Hubble rate, giving rise to a recollapse 
when the singularity is approached. Also a/a and the 
Ricci scalar are bounded and finite. Thus type III singu- 
larities are generically resolved in LQC for all values of 
the curvature index. 

This can be seen from the plots in Fig. [SI where we 
have shown the plots of the Hubble rate and the Ricci 
scalar when the big freeze singularity occurs in the past 
at a — a a . We see that in the classical theory, Hubble rate 
and Ricci scalar diverge as a —> a in the classical theory. 
On the other hand, in LQC the Hubble rate vanishes 
in the past showing that there is a bounce. Further, 
the Ricci scalar remains finite for all the values of the 
curvature index. 

We now make an interesting observation for k = — 1 
case. For a certain choice of parameters which lead to a 
type III singularity in the past, we find that there exists 
an additional classical branch at the small scale factors. 
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FIG. 6: This figure shows the additional branch occurring 
at small scale factors in k — — 1 model for equation of state 
leading to a type III past singularity in Fig. [5] The additional 
branch occurs for certain values of parameters both in the 
classical theory (dashed curve) and LQC (solid curve) when 
scale factor is below the value a , where the big freeze occurs 
classically. Note that, also in this case, only LQC is immune 
from primordial singularity. The parameters used are A — 
— 100, B = -1 and a = 2. 



FIG. 7: Type III future singularity: Classical (dashed) 
and effective LQC (solid) Hubble rate and Ricci scalar with 
k — 0,±1 are compared. In the classical case there is a di- 
vergence of both H and R at a while the curves for LQC 
remain bounded. Notice that LQC cures also the initial sin- 
gularity in the curve case, and in particular the thin (blue) 
solid curves (k — —1) shows a characteristic 'tilt' for small 
a. The parameters are a — 1000, A = 100, B = 1 and a — 2. 



The additional branch faces a big bang singularity in the 
past of classical evolution and is free from singularity in 
future evolution (as is shown by the dashed curves for the 
Hubble rate and the Ricci scalar in Fig. [6]). However the 
additional branch occurs in LQC when the scale factor 
is less than the Planck length. Since the length scale 
involved is below where we expect the effective dynamics 
in LQC to be valid, a more detailed analysis is needed, 
by including modifications pertaining to the inverse scale 
factor, in order to understand the physics emerging from 
LQC in this special case. Never the less, if we assume 
the validity of the effective Hamiltonian in this regime, 
then LQC resolves the past singularity in the additional 
branch (as depicted by the solid curves for the Hubble 
rate and the Ricci scalar in Fig. [6]). 

The type III singularity is also resolved for all values 
of the spatial curvature when it occurs in future. This 
is on the expected lines, as the effects due to intrinsic 
curvature become small at large scale factors. In Fig. [7] 
we depict results from such an evolution. In the classical 
theory, there is a big freeze singularity in future where the 
Hubble rate and the Ricci scalar diverge. For the k = — 1 
universe, there also exists a singularity (big bang) in the 
past evolution. We can see that in LQC there is no type 
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Ill singularity as the Hubble rate vanishes at a = a and 
Ricci scalar reaches a finite value. Further, the past big 
bang singularity in the k = — 1 CclSG, clS pointed above, is 
also resolved. The LQC universe is bounded for all values 
of the curvature index for the equation of state which 
leads to a type III singularity in the classical theory. 

After performing various numerical simulations we 
reach the conclusion that type III singularities, irrespec- 
tive of them occurring in the past or the future are always 
resolved in LQC for all values of the spatial curvature 
index. 




FIG. 8: Type IV past singularity: Hubble rate and Ricci 
scalar for k — 0, ±1. The singularity is at a — 2, where the 
Ricci scalar diverges for k = and k = — 1. Note that for 
k = +1 (dotted line for the classical solutions and thick red 
line for LQC effective solutions) the universe begin at a > a , 
thus the Ricci scalar is always finite. The parameters are 
A = 0.01, B = 1 and a = 1/4. 



4- Type IV singularity: The Big Brake 

Type IV singularity is a derivative-curvature- 
singularity where none of the curvature invariants 
diverge. (In this sense, it does not qualify as a curvature 
singularity). Though the energy density and pressure 
remain finite, a higher derivative of the curvature 
diverges as a — > a Q . The value of a determines the order 
of derivative which blows up [26| . 

As in the type II case, the geodesic equations are well- 
behaved since Hubble rate is finite at a = a . The singu- 
larity is weak because it occurs at a finite volume and the 



theory is geodetically complete since the Hubble rate is 
bounded. Quantum geometric effects have little influence 
on this harmless extremal event beyond which geodesies 
can be extended even in the classical theory [25j . 

The behavior of the curves is very similar to the type 
II case. In Fig. [5] we have shown the results from the 
past singularity. (We have chosen a = 1/4 so that the 
divergence appears in R). We find that the Hubble rate 
vanishes at a = a a both in the classical theory and LQC. 
However, R diverges except for the spatially closed model 
in LQC. This result is similar to what we obtain for type 
II singularity. It turns out that LQC resolves type IV 
weak singularity when occurring in past for the closed 
model for certain choices of parameters. As in the type 
II case, this result highlights the surprises which quanti- 
zation of intrinsic curvature may bring in comparison to 
earlier studies 25]. For k — and k = — 1 model, this 
weak singularity is ignored by the quantum geometry for 
all values of the parameters. 

We show the results for type IV singularity occurring 
in future evolution in Fig.[S] Here we find that for all 
the values of spatial curvature, the Hubble rate vanishes 
at a = a with a divergence in R for both the classical 
theory and LQC. We see that not only the LQC curves 
have no big bang, but this happens also classically for 
the closed universe (Fig.jSJ dotted line). Furthermore, 
for k — +1 it is possible to obtain a baby universe for 
an appropriate choice of the parameters, as for type II 
singularity (see Fig. [3]). In this case R is strongly negative 
but bounded, thus there is no singularity in the past and 
in the future. 



V. SUMMARY 

A fundamental question in quantum gravity is whether 
spacelike singularities of the classical theory are resolved. 
Since not all such singularities signal end of the space- 
time, it is important to understand the role of quantum 
gravitational effects in resolution of strong singularities 
(those beyond geodesies can not be extended) and weak 
singularities (those beyond which geodesies can be ex- 
tended) . These issues were addressed in the loop quanti- 
zation of cosmological models which are spatially flat and 
it was found the non-perturbative loop quantum effects 
resolve all strong singularities and ignore weak singular- 
ities [H. 

The aim of the present analysis was to investigate these 
issues for spatially curved models using phenomenologi- 
cal model of equation of state permitting exotic singulari- 
ties such as big rip, sudden singularities, big freeze singu- 
larity and the big brake singularity. In order to capture 
the role of intrinsic curvature, we considered exotic singu- 
larities both in the future and the past evolution. To our 
knowledge, even in the classical theory exotic singulari- 
ties had not been studies earlier for the spatially curved 
model. For the singularities occurring in the future, the 
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FIG. 9: Type IV future singularity: Comparison of Hub- 
ble rates and R for the classical theory and LQC (in the figure 
LQC curves overlap). The singularity appears in the future 
at a = 1000. The LQC universes (solid lines) do not have 
primordial singularity. Interestingly, for k — 1 this happens 
also classically (dotted line) . In this picture there are no baby 
universe, but they can be obtained for k = 1 for a different 
choice of the parameters. Here the parameters are A = —0.1, 
B — —1 and a = 1/4. 

contribution of the intrinsic curvature is expected to be- 
come very small and we expect results to agree with the 
spatially flat case. This turns out to be true. However, 
more interesting are cases where the exotic singularities 
occur in past. Here one would expect effects due to quan- 
tization of intrinsic curvature to play a non-trivial role. 
In fact we encounter some surprising results. Though 
strong singularities are always resolved in LQC, it turns 



out that for the closed model weak singularities occur- 
ring in the past evolution may also be resolved. Thus 
LQC, does not always ignores weak singularities. This is 
an intriguing result which deserves further investigation. 

Another peculiar feature of the curved models is the 
appearance of a small branch for type II and type IV 
singularities for certain values of the parameters. This 
"baby-universe" is bounded, and is devoid of any singu- 
larities. It will be interesting to analyze these additional 
branches which we find for both spatially open and closed 
models in more details, in particular by taking in to ac- 
count inverse scale factor effects which may play some 
role when scale factor is below the Planck length. 

Our results extend earlier results on generic resolu- 
tion of strong curvature singularities in spatially flat 
model in LQC to the spatially curved models. They also 
bring some important lessons, the primary one being that 
quantization of intrinsic curvature may throw some novel 
unexpected results and we need to gain more insights 
on when quantum gravity may ignore or resolve a weak 
curvature singularity. Another lesson is that as for the 
spatially flat model, spacetime curvature invariants may 
diverge for the spatially curved models and yet there may 
be no physical singularity. 

These results strengthen the case for a generic resolu- 
tion of strong singularities in LQC To achieve this goal, 
the next step will be to include anisotropics and then in- 
homogeneities. The latter will require us to go beyond 
the minisuperspace approximation considered here. This 
brings up additional challenges such as the complete clas- 
sification of the strong and weak singularities in inhomo- 
geneous situations. Two promising avenues where such 
an analysis can be undertaken would be the Gowdy mod- 
els [l6[ and in the spinfoam paradigm ;36] . We hope that 
these studies will also provide insights on the deeper re- 
lation of these frameworks with LQC. 
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